Tunneling in a cavity* 
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Abstract. The mechanism of coherent destruction of tunneling found by Grossmann et al. 
[Phys. Rev. Lett. 67, 516 (1991)] is studied from the viewpoint of quantum optics by considering 
the photon statistics of a single mode cavity field which is strongly coupled to a two-level tunneling 
system (TS). As a function of the interaction time between TS and cavity the photon statistics 
displays the tunneling dynamics. In the semi-classical limit of high photon occupation number n, 
coherent destruction of tunneling is exhibited in a slowing down of an amplitude modulation for 
certain parameter ratios of the field. The phenomenon is explained as arising from interference 
between displaced number states in phase space which survives the large n limit due to identical 



-1/2 



scaling between orbit width and displacement. 



PACS number(s): 42.50.Ct, 42.50.Lc, 05.30.-d 

In recent years the idea of modulation and therefore 
control of tunneling by a monochromatic electromagnetic 
field has been subject of considerable interest. The typi- 
cal Hamiltonian describes a particle in an isolated double- 
well potential (DWP) which is periodically driven by an 
external force 



Hit) = Hy)wp(x) + Sxcos(u>Lt)- 



(1) 



Here TJdwpC^) 1S the Hamiltonian of the DWP, S is the 
amplitude and u>l the driving frequency. The atten- 
tion has mainly focused on a possible enhancement C* 
suppression of coherent tunneling: Lin and BallentineEl 
have demonstrated that the tunneling probability is 
highly enhanced due to periodic modulation for high field 
strengths and driving frequencies close to the classical 
oscillation frequency at the bottom of each well.|-,In the 
opposite limit Grossmann, Hanggi and coworkersa found 
complete suppression of tunneling such that a particle 
initially localized in one of the two wells will never es- 
cape to the other well. They termed this effect "coherent 
destructipn of tunneling" . It has since been of continuing 
interest .ulI The most surprising feature is the periodicity 
of the destruction of tunneling for certain parameter ra- 
tios of S and u>l- So far there is no clear understanding 
of this phenomenon. 

By using the Floquet formalism it has been shown that 
many characteristic features of the tunneling suppression 
can already Jae described in a two-level approximation 
of the DWP.BI Many different aspects of the. effect have 
been illuminated in this framework: MakarovQ, Plata and 
Gomez LlorenteH quantized the electromagnetic field and 
recovered the effect in the limit of iiarge number of pho- 
tons in the field. Wang and ShacQ mapped the driven 
two-level dynamics to a classical one of a charged particle 
moving in a harmonic potential plus a magnetic field in a 
plane. Kayanumafcl explained the suppression of tunnel- 



ing as an effect arising from interference at periodic level 
crossings. 

In the present Rapid Communication, an alternative 
explanation is proposed which uses the concept-af-a/iase 
space interference knowp-fcpm squeezed statesE9iij and 
displaced number statest3'II3 in quantum optics. 

Let us start with a description of the physical situation 
we have in mind. At t — a single two-level tunneling 
system (TS) prepared in a state localized in say the left 
well \L) is injected in a cavity and starts to tunnel be- 
tween its left and right state \L) and \R). The cavity con- 
tains a single mode which has been prepared in a number 
state |n). We consider an ideal cavity, i.e., we neglect any 
kind of dissipation. After an interaction time t{ n the TS 
leaves the cavity and the photon number distribution is 
measured irrespective of the state of the TS. 

Thus our aim will be to calculate the transition prob- 
ability from the product state \L,n) = \L)\n) at t = 
to another product state \i,l) for any i = L, R at time 



t tj 



Pin (ii 



i=L,R 



l(MI*(*m)>| 5 



(2) 



where \^{t in )) = e~ iHU " /h \L,n). For fixed t in this is the 
photon number distribution with P; ra (0) = Si n . For fixed 
I, Pin(Un) describes the dynamics of the cavity mode in- 
teracting with a TS. 

An important aspect of the present problem is the 
strong coupling between the cavity mode and the TS, 
and the separation of timescales between the slow tun- 
neling motion and the fast field oscillations. If we denote 
the coupling energy by g and the tunneling frequency by 
A, this means that we are interested in the limit A <C lol 
and g ~ hear, for low, and \fng ~ frcu^ for high number 
of photons in the field. 
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In this regime the field strongly dresses the TS and 
both TS + field have to be treated as a single unit. In 
contrast to the situation where uil ~ A ^> g/h, we are 
confronted in the present case with a situation in which 
the rotating- wave approximation is not applicable. -Thus 
instead of using the Jaynes-Cummings HamiltonianliJ, we 
must include counter-rotating terms so that our Hamil- 
tonian reads 

H = <j x + ga z (a < + a) + huj^o) a — — (3) 

2 Tilol 

where a', a are the bosonic creation and annihilation 
operators of the cavity mode. We have identified a x — 
\L)(R\ + \R)(L\ and a z = \L)(L\ - \R)(R\. Expressing 
the spin operator in the eigenstates of the TS, |0) = 
2~ 1 /2(| i ) + \R) and |1) = 2- 1 / 2 (|L) - \R), the usual for- 
mulation of the Hamiltionian in quantum optics is recov- 
ered. We also have added a constant energy shift for later 
convenience. It has to be noted that suppression of tun=. 
ncling does not occur in the Jaynes-Cummings model.u 
Hence the situation here is fundamentially different from 
the one usually encountered in quantum optics. 

Owing to the separation of timescales between the tun- 
neling and the oscillation dynamics, PiniUn) can be cal- 
culated in perturbation theory in A/lol by introducing 
dressed states 



|j(n)> 



1 



U\n)-j={\L) + {- 



with j — 0, 1 and 

U = exp[— a z a{o) — a) 



(4) 



(5) 



where a = g/hujL. Equivalently, one may perform the 
polaron transformation H — > IPHU and continue to use 
the product state instead of the dressed state basis. In 
the dressed state basis the Hamiltonian can be written 



as 



H = H D + V 

= \j{m))(mhuj L - (-) J — -^-)0'(™)l 

hA 



^ E b'H)((-) / (i + (-) i " /+m " m ') 

m,7n' 

x £> mm /(2a)(l 

$mm' W(m')\ (6) 



with 



and 



A„ = AD nn (2a) 



(7) 



D ln (a) = (l\D(a)\n) 
1\ 



">e-i\<*\ 2 L l -"(\a\ 2 ) 



where D(a) = exp[a(a^ — a)] is the shift operator and 
L l ~ n (x) an associate Laguerre polynomial (I > n). The 
diagonal part Hd builds a ladder of tunneling dou- 
blets with intra-doublet spacing HA m and inter-doublct 
spacing Hul- Because (l(m)\V\0(m)) = 0, V induces 
only mixing between dressed states belonging to differ- 
ent doublets. Hence, corrections to the dressed states are 
0{A/ujl). Neglecting V for this reason we find for the 
transition probability if / > n 



Pln{U 



£ D lm (a)D mn (-a)e- im " Lti " cos (±A m t, 



D lm (a)D mn (a)e- lM - sin (±A m U 



(9) 



(8) 



and Pi n {tin) — Pni{ti n ) H I < Ti. We conclude that the 
tunneling dynamics can be seen in the spectrum of the 
transition probability of the cavity field. In addition to 
harmonics of ojl its power spectrum also contains reso- 
nances at A rn arising from the tunneling motion. This 
behavior strongly depends on the initial preparation of 
the TS. Injecting it in its ground state |0), for instance, 
yields only resonances at muiL + (— ) ? \ A m where £ = 0,1 
depending on whether n — m is even or odd, respectively, 
and hence only rapid oscillations. 

The effect of a strong coupling between the cav- 
ity mode and the TS is to mix a coherent amplitude 
a = g/TiUL with the intial number state of the mode. 
This happens by displacing the oscillator wave function 
<Pn{x) = (x\n) like x — > x — \^2a in dimensionless coor- 
dinates x — (/iLUL/hy^q. Thus by injecting a TS which 
strongly couples to the caizitjf field, it is possible to realize 
displaced number statesE3l!3. The statistical properties 
of displaced number states have been discussed in Ref. 
p"3| . The photon number distribution is simply given by 
Pr>ns(l) = \Din{o)\ 2 because Di n (a) is the probability 
amplititude of finding I photons in a displaced number 
state \a,n) = D(a)\n). 

The photon distribution ([)]) for number states which 
are dynamically displaced by the tunneling process re- 
sembles Pdns(1)- The displaced number states can either 
be shifted into the same well (first term in (0)) or oppo- 
site wells (second term in (^[)). As expected for a = 0, 
PiniUn) = 5i n independent of the interaction time. 

In Fig. 1(a), (b) we have plotted Pi n {t in ) for n = 10, 
^Ltin = 250 and a = 0.25 and 0.5. Whereas a co- 
herent state_.with n = 10 obeys the familiar Poisson 
distributionEj shown by the line, the dynamically dis- 
placed number states exhibit oscillations shown by the 
histogram. Increasing a results in further oscillations. 
This effect is independent of the specific value of wlUu 
for t in ^ 0, though the absolute value of Pi n (ti n ) depends 
on uiXin- These modulations are clearly exhibited in the 
asymptotic expansion of (^) in the semi-classical limit 
(Bohr's correspondence principle): l,n — > oo, n/l — * 1 
with / — n > finite. If we scale the coupling constant 



2 



between the field and the TS as g oc n -1 / 2 and note that 
associate Laguepe polynomials asymptotically approach 
Bessel functional^ one finds in the semi-classical limit 

Pl-n(tin) = COS 2 (|Atj„) jf_ n Sm.(u> L ti n )\ 

+ sin 2 (|aW) Jf_ n cos(ui L t in ) \ . (10) 

Here 

A = A\J (2Q/lo l )\ (11) 

is the renormalized tunneling frequency, and 

= 2^/ng/h (12) 

is the Rabi frequency. In Fig. 1(c), (d), Eq. (|l"(i| ) is plot- 
ted for H/ujl = 1 and 3.3. It shows the same oscillations 
as the exact expression (||). 

The expression ([To]) is easily understood. The co- 
sine and the sine factors represent the probability for 
the tunneling particle to stay in the well where it has 
been prepared initially, or to escape to the other well, 
respectively.^ The Bessel functions represent the prob- 
ability for the corresponding displaced number state 
to contain I photons if it had n initially (note that 
-Pdns(0 - * Ji- n {^/^L) in the semi-classical limit). 
Hence the Bessel functions represent the photon statistics 
of displaced number states with an effective displacement 
(2Q/u> L ) sin(oj L ti n ) or (2fl/uj L ) cos(u; L t in ). From (0) 
and the cosine and sine factors in (nfl) we further see that 
for specific parameter values of the driving field - where 
20,/ujl hits the roots of the zero order Bessel function 
- the tunneling process is completely suppressed. Note 
that localization cannot occur for small photon number 
in the field because there is no value of a which is si- 
multaneously a root of all Laguerre polynomials. We 
conclude that in the large n limit the tunneling dynam- 
ics is displayed in the photon statistics by an amplitude 
modulation of the cavity mode oscillations. Coherent de- 
struction of tunneling manifests itself in the slowing down 
of this amplitude modulation as depicted in Fig. 2. 

We now show how the present picture can give a quan- 
titative understanding of this effect. First recall that the 
oscillation in the photon number statistics in Fig. 1 orig- 
inate from phase space interference :E3 the nth number 
state can be associated with a circular band of width 
{2n)^ 1 / 2 around its orbit with radius \j2{n + 1/2) and 
centered around the origion. Analogously so can the Zth 
displaced number state which is shifted by \[2ct. Ac- 
cording to the area-of-overlap concepttLJ, the transition 
amplitude between two states is governed by the sum 
of all possible overlap areas weighted with appropriate 
phases. This results in the interference between contri- 
butions from different overlaps which can be constructive 
or destructive, thus, giving rise to oscillations in the pho- 
ton number distribution. Another way of understanding 
the oscillations is that they arise from the possibility of 



a constructive and destructive overlap between displaced 
harmonic oscillator wave functions by noting that Di m (a) 
is a polynomial in I of order min(Z, m), i.e., has min(Z, m) 
roots. 

From this consideration it becomes apparent that vary- 
ing the relative radius of the two bands by changing I vs. 
n, or varying the relative displacement by changing a 
for fixed I and n should result in similar effects. If one 
notices further that the factor Jq{2Vl/ul) in ( |Tl| ) is the 
probability amplitude of finding in a number state dis- 
placed by x — > x — 2^=- again exactly the same number 
of photons, one expects that the periodic suppression of 
tunneling arises from destructive and constructive inter- 
ference of displaced harmonic oscillator wave functions in 
phase space. 

To verify this argument, the n — > oo limit has to be 
considered with care. First take the number of photons in 
the cavity mode to be finite. The factor D nn (2a) which 
scales the doublet splitting (jjj) is the overlap between 
a number state |n) shifted into the right well with the 
same number state shifted into the left well. If we note 
that D nn (2a) is a polynomial in \a\ 2 of degree n, and 
consequently has n zeros, we expect n + 1 oscillations of 
D nn (2a) as a function of a. The n zeros between the 
maxima result from the n possible ways of destructive 
overlap between displaced harmonic oscillator wave func- 
tions. This explains Fig. 2 in Ref. ^. However, contrary 
to the claim there, this picture remains also valid if the 
oscillator energy is larger than the reorganization energy 
(n+ ^)hui l > g 2 /huiL, and even in the limit n — > oo. The 
simple reason for this is that in the semi-classical limit 
n — > oo with 2^/ng/h — ► Q fixed, both the distance be- 
tween the nodes of the harmonic oscillator wave function 
(as a function of x) and the displacement x — > x — 2^=- 

scale as n -1 / 2 . Based on this argument, we expect oscil- 
lation to become noticable as soon as the displacement 
exceeds the bandwidth, i.e., 20,/lol > 1. Finally, we note 
that with the scaling chosen above both displaced orbits 
will always intersect no matter how large SI/lol is. This 
results in the infinite number of oscillations seen in the 
zero order Bessel function ( p"l| ) . 

Summarizing coherent destruction of tunneling is ex- 
plained as a quantum effect arising from the destructive 
interference of displaced harmonic oscillator wave func- 
tions in phase space. The effect is strictly speaking only 
observable in the semi-classical limit for the reasons men- 
tioned above, and survives the large n limit since both 
the bandwidth of the orbit and the displacement scale in 
the same way to zero as n — > oo. 

In closing, we note that experimental study of the ef- 
fects described here is presently still out of range. In a 
microcavity the driving frequency wl is O(GHz) whereas 
typical values of the one-photon Rabi frequency g/Ti are 
0(kHz)Jl3 Hence a number state with large n is needed 
which is difficult to realize experimentally. Possibly a 
trap is more suited because of its lower driving frequency 
lul ~ 0(MHz)Jl3 Damping of the cavity mode which 



3 



is too strong will also render an observation impossible. 
Furthermore the tunneling dynamics must still be coher- 
ent. Finally the superposition state \L) — 2~ 1 / 2 (|0) + 11)) 
is difficult to prepare because the particle experiences 
strong electromagnetic fields when it enters and leaves 
the cavitytZl (for a discussion of this point if one uses 
quantum wells, see Ref. H). 
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FIGURE CAPTIONS 

FIG. 1. Characteristic oscillations in the photon num- 
ber distribution of dynamically displaced number 
states for uilUti = 250, A/lul = 0.2 and |n) = |10). 
In (a) and (b) the quantum limit, Eq. (|^), is dis- 
played for different values of a = g/hujL. In (c) 
and (d) the semi-classical limit, Eq. ([h]), is shown 
for different values of H/ujl = 2^/na. The full line 
shows the Poisson distribution with n = 10. 

FIG. 2. Coherent destruction of tunneling monitored in the 
cavity field. Displayed is Eq. ( |To| ) for I = n, 
A/ujl = 0.2 as a function of LOtMn'- upper figure 
2n/uj L = 2, lower figure 2U/uj l = 2.3. The first 
root of Jq{2VL/ujl) occurs at 2f2/Vk ~ 2.405. The 
effect is exhibited by a decrease of the amplitude 
modulation in the cavity mode oscillations. 
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